Infrared emission in radio galaxy NGC 4261 
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Abstract. We have analyzed the total and nuclear SED for NGC 4261 and find that the dominant 
process for the mid- and far-infrared emission in this object is non-thermal emission from the active 
nucleus. Modeling the emission from the optically detected 300 pc dust disk yields no significant 
disk contribution at any wavelength. To explain the observations, either the disk has an inflated inner 
region which partly absorbs the core, or the intrinsic core spectrum is curved. The inner lOpc of the 
disk can potentially be conceived as an obscuring torus, albeit with optical depth around unity. 



INTRODUCTION 

Unification models for radio-loud active galaxies postulate that for the weaker FR Is only 
the relativistic jet determines the appearance of their nuclei, while in pwerful FR lis a 
circumnuclear dust torus adds to the anisotropy |[lll2]- The infrared emission from this 
torus has been used to test unification models, on the assumption that it dominates the 
infrared output of radio galaxies. 

We have assessed the contribution of the observed nuclear disk to the total infrared 
emission in FRI radio galaxy NGC 4261 (a.k.a. 3C270). This disk is flatter and has 
lower optical depth than a typical obscuring torus, but might be responsible for the 
absence of a distinct UV core. We have used observations from literature and archives, 
paying special attention to the infrared data. When possible, we measured both total and 
unresolved core flux. Between 10 /im and 1 mm NGC 4261 is unresolved. We will use 
nuclear disk for the disk observed in NGC 4261 and torus for the smaller and denser 
torus in FR lis. We adopt Fy °^ v", Hq = 65 km s"^ Mpc^^ and Vsys = 2238 kms. 



ANALYSIS AND MODELING OF THE SED 

In the total flux SED we detect the large scale radio emission and a stellar peak of with 
black body temperature of 4180±40K. Both decline steeply towards the infrared and 
cannot responsible for the far-infrared emission. The core SED shows a typical double 
power-law shape with a turnover around 30 /im, as expected Including only high 
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FIGURE 1. Left : Reference model for NGC 426 1 . The total SED and core SED are shown as asterisks 
and diamonds respectively. Our total model SED is the solid line. The dashed line is the model output 
from dusty disk The dash-dotted line is the disk emission without the stellar heating, and the dotted line 
is the nuclear input spectrum. The model core spectrum is similar to the input spectrum. Error bars are 
omitted where the error is smaller than the symbol. 

Right : The density structure of the disk. Lighter areas represent denser regions. The disk is axis-symmetric 
around the z-axis. 

resolution data, we find spectral indices of 0.2 and -2.5 for the long and short wavelenght 
sides respectively. Including lower resolution VLA data results in a spectral index of - 
0.1, while including the sub-mm data results in 0.1. These values are typical for FRI 
synchrotron cores. 

The 90 jUm detection is well fitted by the power-law, and due to the large calibration 
error is fully consistent with the range of spectral indices we find. Therefore we conclude 
that non-thermal emission is the dominant source for the far-infrared, and find the same 
to be true for the sub-millimeter. The 25 jim detection is consistent with both power-laws 
fit, but due to the large error margin we cannot fully exclude a more significant thermal 
contribution at this wavelength. Our fit results in a sharp and probably unphysical turn- 
over, but fitting a more physical curved spectrum does not alter our findings (see below). 

This implies that the nuclear disk cannot contribute more than 50% of the infrared 
flux. We have tested this with Monte Carlo radiative transfer models, using a setup 
identical to van Bemmel and DuUemond [4] and disk parameters from observations by 
Jaffe et al. |5] and Martel et al. [6]. The dust is heated by both the stellar radiation and 
the active nucleus. We add a UV-X-ray component to the observed core spectrum which 
maximizes the heating from the nucleus. The disk mass, mass distribution and geometry 
are varied within the constraints from observations. From our models we find that the 
disk never contributes more than 10% of the infrared emission (see Fig.[l]). We also find 
that stars dominate the heating of cool dust, while the AGN dominates heating within a 
lOpc radius. However, none of our models provide an explanation for the absence of a 
UV core. Either there is more dust in the line of sight than we account for, or the intrinsic 
core spectrum steepens significantly towards the UV. 
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FIGURE 2. As Fig.[Tl but for the disk with inflated inner region. Here the dash-dotted Hne is the model 
core spectrum. The nuclear input spectrum is flattened in the optical range to improve the fit. 



Within the highly constrained observations, more dust can only be present in the 
innermost region of the nuclear disk, which is unresolved. Evidence for a thicker disk 
within 6pc comes from detections of HI absorption van Langevelde et al. [7]. We 
constructed model disks which have an inflated inner region, with parameters identical 
to the HI disk (see Fig. |2l As expected, these models boost the 2-20 /im emission by 
an order of magnitude, while the far-infrared is unaffected. The inner regions absorb the 
core from the infrared to the UV, but by flattening the intrinsic core spectral index in that 
range, the fit to the observed core spectrum improves and accounts for the lack of a UV 
core. 

Alternatively, if FRIs are mis-aligned BL Lacs, they may have intrinsically curved 
core spectra. Fitting the core spectrum with a model from Chiaberge and Ghisellini [8] 
we can indeed explain the absence of a UV core with a BL Lac type spectrum. All model 
parameters fall well within the observed range for BL Lacs, only the beaming parameter 
is adjusted for the larger inclination. 

Since the observations span about 2 decades in time, core variability can be signifi- 
cant. We therefore allow a factor 2 in our fits, which makes it impossible to distinguish 
between either solution. The 25 and 90 /im detection are also consistent with both, al- 
though there is marginally more room for thermal 25 /im emission with the curved spec- 
trum. However, the curved core spectrum provides much less heating power to the disk 
compared to our models. Therefore, the disk would be less bright, making it difficult 
to account for any thermal emission in the MFIR, and enforcing the conclusion that 
non-thermal emission is the dominant MFIR process. 



CONSEQUENCES FOR UNIFICATION MODELS 



In the FR II unification models, the role of the hypothetical torus is to obscure the nucleus 
at larger inclinations. We have shown that under certain assumptions, the inner parts of 
the nuclear disk are inflated and do exactly that. Classical tori in FRIIs require Ty ^ 1, 
whilst the disk in NGC4261 has Ty ~ 1, so in that sense, it is not an obscuring torus. 
Perhaps we need to adjust our definition and allow tori to have a range of optical depths. 
This would imply that some FR Is have a torus, although it is optically thin, and the inner 
region of the nuclear disk in NGC 4261 would potentially be one. 

If FR Is indeed have thin tori, there may no longer be need for two separate unification 
models. Instead, there may be a continuous distribution of torus optical depths, in which 
FRIs and lis predominantly populate different regions. However, there is not necessarly 
a sharp divide. The SED of the nuclear disk is in good agreement with the infrared 
emission found in broad-line radio galaxies [9]. Assuming their mid-infrared emission 
is thermal, this would imply that some FRIIs could also have thin tori. 

Finally, we would like to stress the importance of using the highest possible resolution 
data to determine the presence of a thermal infrared bump. Lower resolution radio data 
include a significant contribution from older electron populations, which flatten the 
spectrum and thus underpredict the non-thermal infrared emission. In NGC 4261 the 
extrapolation from the total radio spectrum would predict a significant thermal infrared 
bump, which is obviously not the case. In addition, we add to other evidence that the 
far-infrared emission is dominated by stellar heating, which makes it highly unsuitable 
for testing unification models. The best wavelength to detect thermal emission from a 
torus would be between 10 and 20 jim, where both the non-thermal core and the stellar 
emission sharply decline, and the torus emission peaks. 
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